Corridors of barchan dunes: stability and size selection. 
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Barchans are crescentic dunes propagating on a solid ground. They form dune fields in the shape 
of elongated corridors in which the size and spacing between dunes are rather well selected. We 
show that even very realistic models for solitary dunes do not reproduce these corridors. Instead, 
two instabilities take place. First, barchans receive a sand flux at their back proportional to their 
width while the sand escapes only from their horns. Large dunes proportionally capture more than 
they loose sand, while the situation is reversed for small ones: therefore, solitary dunes cannot 
remain in a steady state. Second, the propagation speed of dunes decreases with the size of the 
dune: this leads - through the collision process - to a coarsening of barchan fields. We show that 
these phenomena are not specific to the model, but result from general and robust mechanisms. 
The length scales needed for these instabilities to develop are derived and discussed. They turn out 
to be much smaller than the dune field length. As a conclusion, there should exist further - yet 
unknown - mechanisms regulating and selecting the size of dunes. 

PACS numbers: 45.70.Qj; 45.70. Vn; 89.20.-a 



Since the pioneering work of Bagnold [lj , sand dunes 
have become an object of research for physicists. Ba- 
sically, the morphogenesis and the dynamics of dunes 
result from the interaction between the wind, which 
transports sand grains and thus modifies the shape of 
the dune, and the dune topography which in turn con- 
trols the air flow. A lot of works have been devoted to 
the study of the mechanisms at the scale of the grain 
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The interested reader should refer to a previous paper 
p2| for a review of these works. Our aim is to focus here 
on dune fields and to show that most of the problems at 
this scale are still open or even ill-posed. 

The most documented type of dune, the barchan |46| . 
is a crescentic shaped dune, horns downwind, propagat- 
ing on a solid ground. In the general picture emerg- 
ing from the literature, barchans are thought as soli- 
tary waves propagating downwind without changing their 
shape and weakly coupled to their neighborhood. For 
instance, most of the field observations concern geomet- 
ric properties (morphologic relationships) and kinematic 
properties (propagation speed). This essentially static 
description probably results from the fact that barchans 
do not change a lot at the timescale of one field mission. 

As shown on figure 2] barchans usually do not live iso- 
lated but belong to rather large fields [IJj ■ Even though 
they do not form a regular pattern, it is obvious that 
the average spacing is a few times their size, and that 
they form long corridors of quite uniformly sized dunes. 
Observing the right part of figure ^ the barchans have 
almost all the same size (6 m to 12 m high, 60 m to 
120 m long and wide). Observing now the left part of 
figure n the barchans are all much smaller (1.5 m to 3 m 





FIG. 1: Aerial photograph showing part of the barchan field 
extending between Tarfaya, Laayoune and Sidi Aghfinir in 
southern Morocco, former Spanish Sahara. The trade winds, 
dominant in the region, blow from the north (from the top 
of the photograph). Several corridors are visible in which 
the size of barchans and their density is almost uniform. As 
confirmed by the three zooms, the size of dunes is different 
from one corridor to another. 
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high, 15 m to 30 m long and wide) and a small band can 
be distinguished, in which the density of dunes becomes 
very small. Globally, five corridors stretched in the di- 
rection of the dominant wind can be distinguished: from 
right to left, no dune, large dunes, small dunes, no dune 
and small dunes again. Figure^shows only 17 km of the 
barchan field. Direct observations show that these five 
corridors persist in a coherent manner over hundreds of 
kilometers along the dominant wind direction. 

The content of this paper is perhaps a bit unusual as 
we will mostly present negative results. Indeed, we will 
show that none of the dune models [3^l3^.l35ll3a.[ 3 ^, l33. 
[MOD. 03 nor the coarse grained field simulation [44| are 
able at present to reproduce satisfactorily the selection 
of size and the formation of corridors. 

More precisely, we shall first address the stability of 
solitary dunes, and conclude that, given reasonable or- 
ders of magnitude for dune sizes and velocities, barchans 
which are considered as "marginally unstable" by other 
authors |4l| would in fact have the time to develop their 
instability over a length much smaller than that of the 
corridor they belong to. As a consequence, isolated dunes 
must be considered as truly unstable objects. Further- 
more, the origin of this instability is rather general and 
model independent, as it can be understood from the 
analysis of the output sand flux as a function of the dune 
size. 

One can wonder whether interactions via collisions be- 
tween dunes can modify the dynamics and the stability 
of dunes. In a recent paper Lima et al. [44j have in- 
vestigated the dynamics of a field and have claimed to 
get realistic barchan corridors. However, they made use 
of numerical simulations into which individual dunes are 
stable objects of almost equal size (6% of polydispersity) . 
They consequently obtained a nearly homogeneous field 
composed of dunes whose width is that of those injected 
at the upwind boundary. We show here that the actual 
case of individually unstable dunes leads by contrast to 
an efficient coarsening of the barchan field. 

The paper is organized as follows. In order to get a 
good idea of the mechanisms leading to these two in- 
stabilities, we first derive a 3D generalization of the C£ 
model previously used to study 2D dunes [3] . We then 
show that the two instabilities predicted by the C£ model 
are in fact very general and we will derive in a more gen- 
eral framework the time and length scales over which they 
develop. Turning to field observations, we will conclude 
that the formation of nearly uniform barchan corridors is 
an open problem: there should exist further mechanisms, 
not presently known and may be related to more compli- 
cated and unsteady effects such as storms or change of 
wind direction, to regulate the dune size. 



I. BARCHAN MODELING. THE Q MODEL 

We start here with the state of the art concerning the 
modeling of dunes by Saint- Venant like equations. First, 



the mechanisms of transport at the scale of the grain H0, 
II II II II II 13 IH 13 Hi El 13 El 13112 determine 
at the macroscopic scale - at the scale of the dune - 
the maximum quantity of sand that a wind of a given 
strength can transport. As a matter of fact, when the 
wind blows over a flat sand bed, the sand flux increases 
and saturates to its maximum value Q after a typical 
length L called the saturation length 1, 16, 28, 29]. This 
length determines the size of the smallest propagative 
dune. 

The other part of the problem is to compute the tur- 
bulent flow around a huge sand pile of arbitrary shape 
[30l l3ll | . Since the Navier-Stokes equations are far too 
complicated to be completely solved, people have derived 
simplified descri ptio ns of the turbulent boundary layer 

pi in hi m m m m m eu. The first s t eP 

initiated by Jackson, Hunt et al. has been to derive an 
explicit expression of the basal shear stress in the limit of 
a very flat hill. Kroy et al. [3 EH 

have shown that this 
expression can be simplified without loosing any impor- 
tant physical effect. In particular, it keeps the non-local 
feature of the velocity field: the wind speed at a given 
place depends on the whole shape of the dune. 

Being a linear expansion, this approach can not ac- 
count for boundary layer separation and in particular for 
the recirculation bubble that occurs behind dunes. Fol- 
lowing Zeman and Jensen |3l| and later Kroy et al, the 
Jackson and Hunt formula is in fact applied to an enve- 
lope of the dune constituted by the dune profile prolonged 
by the separation surface. 

As already stated in one of our previous papers [3] , we 
proposed to name C£ the class of models which describe 
the dynamics of dunes in terms of the dune profile h and 
the sand flux q, and which include (i) the mass conser- 
vation, (ii) the progressive saturation of sand transport 
and (Hi) the feedback of the topography on the sand ero- 
sion/deposition processes. We chose this fancy name in 
reference to the spatial organization of the dunes which 
propagate like the flight of wild ducks and geese. 



A. 2D and 3D main equations 



Let us start with a quick recall of the set of 2D C c 
equations that we already introduced in |28|. Let x de- 
note the axis oriented along the wind direction, and t 
the time. The continuity equation which ensures mass 
conservation simply reads 



d t h + d x q = 0. 



(1) 



Note that q(x, i) denotes the integrated volumic sand 
flux, i.e. the volume of sand that crosses at time t the po- 
sition x per unit time. The saturation process is modeled 
by the following charge equation 



d x q 



Qsat 



L 



(2) 
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It is enough to incorporate the fact that the sand flux 
follows the saturated flux q sat with a spatial lag L. It is 
a linearized version of the charge equation proposed by 
Sauermann et al. fl6j ]. 

The saturated flux q sat is a growing function of the 
shear stress. This shear stress can be related to the dune 
profile h by the modified Jackson and Hunt expression. 
Since this expression comes from a linear expansion, we 
can directly relate q sa t to h by: 

S^M = 1 + A f^d x h e (x- X ) + Bd x h e (x), (3) 

Q J n x 

where Q is the saturated flux on a flat bed and h e the 
envelope prolonging the dune on the lee side (see Ap- 
pendix and |H, [5j for the details of construction) . The 
last term takes into account slope effects, while the con- 
volution term encodes global curvature ones. The only 
relevant length scale is the saturation length L of the 
sand flux. The other relevant physical parameter is the 
saturated sand flux on a flat bed Q. All the lengths are 
calculated in units of L, time in units of Q/ L 2 , and fluxes 
in unit of Q. A and B could in principle be predicted by 
the Jackson and Hunt analysis but we rather take them 
as two tunable phenomenological constants. 

In three dimensions, equations are very similar, albeit 
slightly different. In order to express the total sand flux 
(which is now a 2D vector), we need to distinguish saltons 
and reptons The reason is that in contrast to the 

saltons, which follow the wind, the motion of the reptons 
is sensitive to the local slope 3] . Because the reptons are 
dislodged by the saltons, we assume that their fluxes are 
proportional [Tt| . so that the total flux can be written 
as the sum of two terms, one along the wind direction x 
and the other along the steepest slope 0: 

qtot = qx- DqVh. (4) 

The continuity equation then takes its generalized form 

d t h + V-q tot =0. (5) 

The down slope flux of reptons acts as a diffusive process. 
The diffusion coefficient is proportional to q so that no 
new scale is introduced - D is a dimcnsionlcss parameter. 
This diffusion term introduces a non-linearity that has a 
slight effect only: almost the same dynamics is obtained 
if a constant diffusion coefficient is used instead. Equa- 
tions J2| and © can be solved independently in each 
slice along x. 

In summary, the C£ model considered here includes in 
a simple way all the known dynamical mechanisms for 
interactions between the dune shape, the wind and the 
sand transport. 

B. Propagative solutions of the Cc model. 

In a previous paper [2^, we have studied in details 
steady propagative solutions in the 2D case. They also 




FIG. 2: Evolution from a conical sand pile to a steady prop- 
agative barchan computed from the Cc model. To obtain this 
steady solution the output flux is re-injected homogeneously 
at the upwind boundary. Times are given in units of L 2 /Q. 
Stereoscopic view: a) place the figure at ~ 60 cm from your 
eyes b) focus behind the sheet, at infinity (you should see 
three dunes) c) focus on the middle dune and relax d) you 
should see the shape in 3D. 
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FIG. 3: Relationship between the inverse velocity 1/x and the 
width w for barchans in the steady state. The line corresponds 
to the best fit by a Bagnold-like relation of the form x = 
aQ/(w + w c ). It gives a — 56 and w c = 9.5. 
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FIG. 4: Relationship between the volume V and the width w 
of solitary barchans. Note the log-log scales. The solid line 
corresponds to the best fit by the relation V = bw 2 (w + w v ). 
It gives b = 0.011 and w v = 22.9. 



apply to transverse dunes, i.e. invariant in the y direc- 
tion. We will now focus on three dimensional solitary 
dunes computed with the model presented in the lat- 
ter section. The details of the integration algorithm and 
the numerical choice of the different parameters can be 
found in the Appendix and a more detailed discussion 
about the influence of the diffusion parameter is discussed 
in |45| . Figure |21 shows in stereoscopic views the time 
evolution of an initial conical sand pile (t = 0). Horns 
quickly develop (t — 16 and t — 32) and a steady barchan 
shape is reached after typically t — 50. Note that the 
propagation of the dune is not shown on figure [21 the 
center of mass of the dune is always kept at the center of 
the computation box. 

The ori gina l Cq model proposed by Kroy, Sauermann 
et al. |40l Elj was the first of a long series of models in 
which a steady solitary solution could be exhibited, with 
all the few known properties of barchans. In particular, 
the dunes present a nice crescentic shape with a length, 
a width, a height and a horn size that are related to each 
others by linear relationships. They propagate downwind 
with a velocity inversely proportional to their size, as 
observed on the field. These properties are robust inside 
the class of modeling, since we get the same results with 
the simplified version that we use here. We will only 
show in the following two of these properties, important 
for the stability discussion, namely the velocity and the 
volume as functions of the dune size. 

Since they are linearly related one to the others, all 
the dimensions are equivalent to parameterize the dune 
size. We choose the width w as it is directly involved in 



the expression of the sand flux at the rear of the dune. 
Figure |3 shows the inverse of the propagation velocity of 
the dune as a function of w. The velocity decreases as 
the inverse of the size: 



Preliminary fields measurements of the displacement of 
the dunes shown on figure ^ over 27 years have given 
aQ = 3700 m 2 /year and w c = 33 m. The transverse 
velocity y is found to be null, as lateral inhomogeneities 
of the sand flux are unable to move dunes sideways |44| . 
The volume of V is plotted on figure 0] This relation is 
well fitted by: 

V = bw 2 (w + w v ), (7) 

where the numerical coefficients are b ~ 0.011 and w v ~ 
22.9. This value roughly corresponds to the volume of 
a half pyramid, with a height h ~ 0.1 w and a width w 
which gives a volume V ~ w 3, /60. One can observe that 
barchan dune are not self similar object: the deviation 
observed for small dunes is related to the change of shape 
due to the existence of a characteristic length L. 



C. Instabilities 

The choice of the boundary conditions is absolutely 
crucial: to get stationary solutions, the sand escaping 
from the dune and reaching the downwind boundary is 
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uniformly re-injected at the upwind one. Obviously, this 
ensures the overall mass conservation. Doing so, the sim- 
ulation converges to a barchan of well defined shape of 
width Woo with a corresponding sand flux q^. 

However, under natural conditions, the input flux q is 
imposed by the upwind dunes. We thus also performed 
simulations with a given and constant incoming flux. Fig- 
ure [5] shows the evolution of two dunes of different sizes 
under an imposed constant input flux. One is a bit larger 
than the steady dune corresponding to the imposed flux, 
and the other is slightly smaller. It can be observed that 
none of these two initial conditions lead to a steady prop- 
agative dune: the small one shrinks and eventually dis- 
appears while the big one grows for ever. The steady 
solution obtained with the re-injection of the output flux 
is therefore unstable. 

If solitary dunes are unstable, it is still possible that the 
interaction between dunes could stabilize the whole field. 
It is not what happened in the Cq model. Instead, an 
efficient coarsening takes place as shown by Sauermann 
in the chapter 8 of |4l|. 

As a first conclusion, the Cq model predicts that soli- 
tary barchans and barchan fields are unstable in the case 
of a permanent wind. We will see below that these two 
instabilities are generic and not due to some particularity 
of the modeling. In particular, they can be also observed 
with the more complicated equations of Kroy, Sauermann 
et al. who deal with a non-linear charge equation and 
take explicitly into account the existence of a shear stress 
threshold to get erosion. 

Therefore, we can wonder what are the dynamical 
mechanisms responsible for these instabilities. Would 
they have time/length to develop in an actual barchan 
field? Seeking answers to these questions, we will now 
investigate the two instabilities in a more general frame- 
work. As a first step, we will investigate the time 
and length scales associated to the evolution of barchan 
dunes. 



the turnover time: 



In the geological community, the turnover time is be- 
lieved to be the time after which the dune looses the 
memory of its shape. The idea is that a grain remains 
static inside the dune during a cycle of typical time Tt : it 
then reappears at the surface and is dragged by the wind 
to the other side of the dune. In other words, after Tt all 
the grains composing a dune have moved, and the inter- 
nal structure of the dune has been renewed. But this does 
not preclude memory of the dune shape at times larger 
than Tt, and one can wonder whether Tt is the internal 
relaxation time scale to reach its equilibrium shape. The 
scaling JBJ of the propagation speed involves the cut-off 
length scale w c , which can be measured by extrapolat- 
ing the curve of|3|to zero. Note that the existence of a 
characteristic length scale also appears in the dune mor- 
phology |2£| |4(J, El[ • In the following, we will assume 
that the barchans are sufficiently large to be considered 
in the asymptotic regime. We checked that introducing 
cuts-off w c or w v to capture the shape of curves like that 
of figures |3 or 0] in the region of small w does not change 
qualitatively the results. In the following we then take 
w c = and w v — for simplicity. Under this assump- 
tion, using the expression © of the propagation speed, 
the turnover time reads: 

2 

w 

Tt = -Ty (9) 

aQ 

Of course, the length scale At associated with the 
turnover time is the size of the dune itself: 

A t = w. (10) 
B. Relaxation time 



II. TIME AND LENGTH SCALES 

Three different time scales govern the dynamics of 
dunes: a very short one for aerodynamic processes (i.e. 
the grain transport), the turnover time for the dune mo- 
tion, and a much larger time scale involved in the evolu- 
tion of the dune volume and shape under small pertur- 
bations of the wind properties. 



A. Turnover time 

The dune memory time is usually defined as the time 
needed to propagate over its own length. Since the length 
and the width of the dune are almost equal - this is only 
a good approximation for steady dunes - we will use here 



Let us consider, now, a single barchan dune submitted 
to a uniform sand flux. The evolution of its volume is 
governed by the balance of incoming (j) in and escaping 
4>out sand volumes per unit time: 

V = <p in - (j)out- (11) 

4>i„ is directly related to the local flux q upwind the 
dune, defined as the volume of sand that crosses a hor- 
izontal unit length line along the transverse direction y 
per unit time. Assuming that this flux q is homogeneous, 
the dune receives an amount of sand simply proportional 
to its width w: 

(pin = qw. (12) 

The loss of sand 4> Q ut is not simply proportional to 10 
because the output flux is not homogeneous. Figure 
shows the flux in a cross-section immediately behind the 
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FIG. 5: Origin of the flux instability. Two dunes are submitted to a given flux, which is not the equilibrium sand flux for 
both dunes. The small dune (dotted line) is then under supplied an can only shrink. On the contrary, the bigger dune (solid 
line) receive too much sand and grows. The evolution of their volume V is depicted on the right. Eventually, the small one 
disappears. 
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FIG. 6: Top: three dimensional shape of a barchan dune 
obtained with the C£ model, with semi-periodic boundary 
conditions to ensure the mass conservation. Bottom: profile 
of the resulting output sand flux. The sand loss is localized 
at the tips of the horns. There, the flux is almost saturated: 



FIG. 7: Output flux 4> ou t as a function of the barchan width 
w for the equilibrium input flux (white circles) and for a null 
input flux (black circles). <f> ut is not simply proportional to 
w and does not vanish at small size. As a consequence, the 
sand loss is proportionally smaller for a large dune than for 
a small one. The straight line corresponds to the best fit by 
a function of the form (j> ut = Q(A + aw). It gives A = 1.33 
and q = 0.05. 



dune. One can see that the sand escapes only from the 
tip of the horns, where there is no more avalanche slip 
face. As a matter of fact, the recirculation induced be- 
hind the slip face traps all the sand blowing over the 
crest. We computed in the model the output flux 4> out 
as a function of the dune width w (figure [TJ. Within a 
good approximation it grows linearly with w: 

<l>out ~ Q(A + aw). (13) 

For the set of parameters chosen, the best fit gives 



a = 0.05 and A = 1.33 L. Note that the discrepancy 
°f 4>mtt with the linear variation for small dunes can be 
understood by the progressive disappearance of the slip 
face (domes), leading to a massive loss of sand. 

It can be observed from figure El that q is almost sat- 
urated in the horns. The ratio 4> ou t/Q then has a geo- 
metrical interpretation as it gives an estimate of the size 
of the horn tips. Therefore, in the C£ model, the horn 
size is not proportional to the dune width, but grows as 
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A + aw. This is consistent with the observations made by 
Sauermann et al. in southern Morocco: they claim that, 
at least for symmetric solitary dunes, the slip face is pro- 
portionally larger for large dunes than for small ones, i.e. 
that the ratio of the horns width to the barchan width 
decreases with w. 

With these two expressions for the input and output 
volume rates, the volume balance reads: 



V = qw - Q(aw + A). 



(14) 



If we call Woe and goo the width and the flux of 
the steady dune for which the dune volume is constant 
(V = 0), we can define r r — (V — Voo)/V, taken around 
the fixed point. We get: 



3^4 
QA 



(15) 



It also gives us the relaxation length for the dune A r , 
which is the distance covered by the dune during the time 



i.e: 



\ r — 



Zabw 2 . 



(16) 



C. Flux screening length 

For a dune field, the situation is a little bit more com- 
plex. The flux at the back of one dune is due to the 
output flux of an upwind dune. The latter is strongly 
inhomogeneous since the sand is only lost by the horns 
(figure HJJ ■ Field observations show that there is a sand 
less area downwind of the barchans - see also the inset 
of figure |H| This zone is larger than the recirculation 
bubble and indicates a small amount of sand trapped by 
the roughness of the ground. The fact that this 'shadow' 
heals up is a signature of a lateral diffusion of the sand 
flux. The length of the shadow is typically a few times 
the dune width and is in general smaller than the dis- 
tance between dunes. So, the flux can be considered as 
homogeneous when arriving at the back of the next dune. 

The distance X q over which the flux changes is thus the 
distance, along the wind direction, between two dunes. 
It is the mean free path of one grain traveling in straight 
line along the wind direction. Let us consider an homo- 
geneous dune field composed of identical dunes of width 
Woo. The number of dunes per unit surface is Noo- It 
can be inferred from figure[H]that on average there is one 
dune in the surface X q Woo (colored in gray on figure |SJ. 
The flux screening length thus depends on the density of 
dunes as: 



1 



NocWo 



(17) 



Note that this length is larger than the average dis- 
tance (iVoo) -1 / 2 between dunes - just like the mean free 
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FIG. 8: The flux screening length is the mean free path 
along the wind direction. In other words, it is the mean lon- 
gitudinal distance between two dunes. Inset: the sand flux is 
much larger on the back of the dune (dark zone) than on the 
surrounding ground (gray zone). Downwind the dune, it be- 
comes inhomogeneous and in particular, it is null inside the 
recirculation bubble (white zone) and low in the triangular 
shadow zone (light grey). After few dune sizes, the diffusion 
of grains rehomogenizes the flux. 



path in a gas. Since the grains in saltation on the solid 
ground go much faster than the dune (by more than five 
orders of magnitude), the flux screening time r q can be 
taken as null: 



0. 



D. Orders of magnitude 



(18) 



These different time scales can be estimated using the 
orders of magnitude obtained from field observations in 
the region of figure ^ The velocity/width relationship 
has allowed to estimate aQ ~ 3700 m 2 /year and w c ~ 
33 m in this region. Combined with C£ model results, this 
gives estimates of the saturation length L — 3.5m, the 
minimal horn width A = 4.6 m and the saturated flux 
far from any dune Q ~ 66 m 2 /year. These values are 
corroborated by direct measurements of these quantities 

Let us consider a small dune of width 20 m and a 
large dune of width 100 m belonging to the corridors of 
dunes shown on figure ^ The distance \ r covered by 
the dune before the equilibrium between the size and the 
sand flux be reached is respectively 160 m and 4 km. 
In all the cases, it is much smaller than the dune field 
extension (typically 100 km corresponding to 5000 small 
dune widths or 1000 large dune widths). Obviously, A,. 
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is much larger than the turnover length At = w, and it is 
therefore clear that the turnover scales does not represent 
the memory of the dune. 

The density of dunes can be inferred from figure ^ arid 
is around 0.1 /w^ (the average distance between dunes 
is around 3 dune sizes) . Directly from figure ^ or from 
formula (|17l) . the flux adaptation length X q is around 10 
dune sizes i.e. 200 m for the small dune and 1 km for 
the large one. Obviously, X q can be very different from 
place to place. For instance, the left corridor shown on 
figure n is much denser than the third from the left. If 
the density of small dunes is 1 /w^ instead of 0.1 /w^, 
X q becomes equal to the dune size (20 to). 

Using the previous value of Q, the dune velocities are 
180 m/year and 37 m/year for the 20 m and 100 m 
barchans respectively. The corresponding turnover times 
Tt are 5.2 weeks and 2.7 years, while the relaxation time 
r r is as large as 10 months for the small dunes and 
1.1 century for the large ones. Finally, the flux adap- 
tation time r q is equal to the flux screening length A,j 
divided by the grain speed (~ 1 m/s). It can thus be es- 
timated to 3 minutes for 20 m barchans and 16 minutes 
for the 100 to ones. 

The scale separation of the three times is impressive. 
r q ~ 3 minutes <C r t ~ 5.2 weeks <C r r ~ 10.4 months 
for the small dunes, whereas for the large ones it reads 
r q ~ 16 minutes <C r t ~ 2.7 years <C r r ~ 1.1 century. 
This shows that the annual meteorological fluctuations 
(wind, humidity) have potentially important effects: the 
actual memory time is always larger than seasonal time. 

Sauermann has estimated a characteristic time for the 
evolution of the volume of a 100 m wide dune. He found 
several decades |41j , which is comparable to the value we 
found for r r . On this basis he concluded that " consider- 
ing this timescale it is justified to claim that barchans in 
a dune field are only marginally unstable." This is a mis- 
leading conclusion as the length A r (4 km for w = 100 m) 
should be compared to the corridor size (~ 100 km). 
Moreover, for small dunes as those on the left corridor of 
figure^ A r is found to be as small as 160 to which is of the 
order of one hundredth of the portion of field displayed 
on the photograph. The evolution time and length scales 
could be thought as 'very large' - with respect to human 
scales - but compared to the dune field size, they turn 
out to be small. Therefore, barchans have in fact the 
time and space to change their shape and volume along 
the corridors. 



III. FLUX INSTABILITY 

A. Stability of a solitary barchan 

We seek to understand the generality of the instabili- 
ties revealed by the C£ model. We will first investigate 
theoretically the stability of a solitary dune in a constant 
sand flux q. We recall that the overall volume of sand re- 
ceived by this dune per unit time is simply proportional 
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FIG. 9: The output volume rate (pout = Q(A + aw) gives, 
when compared to the input rate pin = <l w , steady solutions 
that are unstable if A > 0, and stable if A < 0. The main 
point is to see whether the two lines cross from below or above 
at the steady point. 



to its width: (p in — qw. As found in the C c model, we 
suppose that <p out = Q(A + aw) (figure left). Let us 
investigate what happens for different values of the input 
flux q. 

If q < a (dot-dashed line) the two curves (pi n {w) and 
<t>out{w) do not cross, which means that no steady so- 
lution can be found. Since the input sand volume rate 
is too low, any dune will shrink and eventually disap- 
pear. On the other hand, a fixed point does exist 
for q > a (thin solid line). Suppose that this dune is 
now submitted to a slightly larger (resp. smaller) flux q 
(dotted lines): it will grow (resp. shrink). However, the 
corresponding steady states are respectively smaller and 
larger, so that they cannot be reached dynamically. We 
now fix the input flux to g^, and change the dune size w, 
as in figure A dune of width slightly smaller than w^ 
under will shrink more and more because it looses sand 
more than it earns. In a similar way, a dune larger than 
Woo will ever grow. In other words, the steady solutions 
are unstable. 

This mechanism explains the flux instability of Cq 
barchans. This stability analysis is in fact robust and 
not specific to the linear choice for (pout- Any more com- 
plicated function would lead to the same conclusion pro- 
vided that (pi n crosses (pout from below. The stability only 
depends on the behavior of the <p's in the neighborhood 
of the steady state. 

How could a solitary barchan be stable? It is enough 
that (pi n crosses (p out from above. Without loss of gen- 
erality, we can keep a linear dependence of <p ou t on w in 
the vicinity of the fixed point, but this time with A < 
(figure right). In this case, the situation for which the 
input sand flux q is larger than a (dot-dashed line) leads 
to an ever growing dune. Steady solutions exist when 
q < a. Because a smaller (resp. larger) sand flux now 
corresponds to a smaller (resp. larger) dune width, these 
solutions are, by contrast, stable. 
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FIG. 10: Sketch showing the instability due to the exchange 
of mass between the dunes: dunes slightly smaller than Woo 
loose sand and make their downwind neighbors grow. Note 
that the dune field is assumed to remain locally homogenous. 



In a more quantitative and formal way, the mass bal- 
ance for a barchan (|14(l can be rewritten in terms of the 
dune width only: 



Q(A + aw) 



(19) 



we 



3bw 2 

Linearizing this equation around the fixed point w c 
obtain: 

T r W = W — Woo. (20) 

The sign of the relaxation time t t is that of A - see 
relation (|15|) . Therefore if A is positive, w will quickly 
depart from its steady value Woo. In the inverse case 
A < 0, any deviation of w will be brought back to Wk,. 

In summary, the stability of a solitary barchan depends 
whether the ratio of the output volume rate to its width 
4>out/w increases or decreases with w. This quantity is 
perhaps not easy to measure on the field but we have 
shown that it is directly related to the ratio of the size of 
the horn tips to the dune width. If viewed from the face, 
the horn tips become in proportion smaller as the dune 
size increases, the barchan is unstable. This is what is 
predicted by the G~ model, in agreement with the few 
field observations 



B. Stability of a dune field 

At this point, the stability analysis leads to the fact 
that a single solitary barchan is unstable. Could then 
dunes be stabilized by their interaction via the sand flux? 
Let us consider a dune field which is locally homogeneous 
and composed, around the position (x,y), barchans of 
width w with a density N. We ignore for the moment 
the fact that dunes can collide. The conservation of the 
number of dunes then reads 



N = d t N + d x (±N) + d y (yN) = 0. 



(21) 



FIG. 11: Sketch showing the instability due to the collisions 
between the dunes. If one dune is slightly larger than the 
others, it goes slower and will absorb incoming dunes. 



The eulerian evolution of the dunes width is given by: 

w = d t w + xd x w + ijdyW. (22) 

This equation has to be complemented by the equation 
governing the evolution of the sand flux q between the 
dunes which results from the variation of the volume of 
these dunes (and reciprocally): 



8 x q = -NV. 



(23) 



These three equations can be closed using the previ- 
ous modeling of x, y = 0, w and V = 3bw 2 w. Any 
homogeneous field of barchans of width and den- 
sity iVoo is a solution, provided that there is a free flux 
<Zoo = Q(A/uioo + a) between the dunes. 

We are now interested in the stability of this solution 
towards locally homogeneous disturbances. Note that 
such a choice is still consistent with equations that do 
not take collisions into account. We thus expand w,q 
and N around their stationary values , , , and 
introduce the length and time scales A r ,A,j and r r . We 
get: 



r r d t N + X r d x N — ^^\ r d x w = 0, 

Woo 



(24) 



r r d t w + X r d x w = (w -Woo) + — ^(q - q^), (25) 

\d x q = -{q - qoo) - ^-(w ~ Woo)- (26) 

w 

Without loss of generality we can write the disturbances 



n „<rt+ikx 



, w-w c 



at+ikx 



under the forms: q—qco — QiS 
and iV — iV^ = Nie at+lkx . Solving the system of linear- 
equations we obtain the expression of the growth rate a 
as a function of the wavenumber k: 



= ik 



i + (k\ q y 



-A r + 



(fcA,) 2 



(27) 



The sign of the real part of the growth rate a is that of r r 
and thus of A. The stability of the dune field is therefore 
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that of the solitary dune. If A < 0, which means that all 
the individual dunes are stable, the field is (for obvious 
reasons) stable. But in fact, all the individual dunes are 
unstable ( A > 0) , so that a field in which dunes interact 
via the sand flux is also unstable. 

The result of the above formal demonstration, can also 
be understood via a simple argument illustrated on fig- 
ure llOl Consider a barchan dune field at equilibrium: for 
each dune, input and output volume rate are equal. Now, 
imagine that the input flux of a dune slightly decreases 
for some reasons. As explained in the previous subsec- 
tion, if this dune is unstable (A > 0) it tends to shrink. 
Consequently its output flux increases, and makes its 
downwind neighbors grow. Therefore, even a small per- 
turbation of the sand flux can dramatically change the 
structure of the field downwind. 



IV. COLLISIONAL INSTABILITY 
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The free flux is not the only way barchans can influ- 
ence one another. If sufficiently close, they can interact 
through the wind, i.e. aerodynamically. This is possi- 
ble when the dunes get close to each other. In this case, 
they actually collide. We therefore would like to investi- 
gate the behavior of one particular dune in the middle of 
the field. 

Let us consider a homogeneous field of barchans of 
width Woo, with an additional dune of size w = (l+?7)ifoo- 
The variation of the volume of this dune is due to the 
sand flux as well as the collisions of incoming dunes. 
These collisions are a direct consequence of the fact that 
smaller dunes travel faster (Eq. [fjj). The number of col- 
lisions per unit time is proportional to the dune density 
Noo times the collisional cross section w + Woo times the 
relative velocity aQ (l/w^ — 1/w). We assume that the 
collisions lead to a merging of the two dunes. Then, each 
collision leads to an increase of the mass of the larger 
dune by I4o = bw^. We can then write for this particu- 
lar dune: 

V = q W - Q(aw + A) 

„ aQ(w — Woe) 

UTOoo 

Introducing a critical dune density N c as 

-A 



N r = 



2abw3 ' 



(29) 



the equation governing the evolution of the width of the 
dune considered reads: 



V 



(2 + ??) N a 
2(1 + 7?) N c ) (1 + tj) 2 ' 



(30) 



Note that, rigorously speaking, N c is a positive quantity 
and thus a true density only for A < (see below). Fig- 
ure 1 121 shows rj as a function of the dune size. Expanding 



FIG. 12: Growth rate r r r\ of a dune due to collisions, as a 
function of the rescaled size rj in three different cases. If the 
solitary dune is unstable (A > 0), the field is also unstable 
towards the collisional instability (dot-dashed line). If the 
solitary dune is stable (A < 0), the stability of the field de- 
pends on the dune density. At high density (solid line), the 
field is linearly unstable while at low density (dashed line) , it 
is stable towards any disturbance. 



linearly around r\ = 0, we obtain the growth rate a = fj/rj 
as: 



T r (J = 1 



^Voc 
N r . 



(31) 



Therefore, if the dunes are individually unstable (case 
A > 0, figure El a )i the dunes are always unstable 
towards the collisional instability. The barchan field 
quickly merge into one big barchan dune. If the dune are 
individually stable (A < 0), the same instability develops 
but only when the dune density is larger than the critical 
dune density N c (figure IT21 bl . Suppose indeed that one 
collision occurs in the middle of an homogeneous field, 
creating a dune of twice its original volume. Since it is 
larger, this dune slows down and a second collision oc- 
curs before the large dune has recovered its equilibrium. 
If now the dune density is small (figure El c ) > the time 
before a second collision happens is sufficiently large to 
allow the large dune to recover its equilibrium volume, 
and in this case a field of stable dunes is stable towards 
the collision process. 



V. BARCHANS CORRIDORS, AN OPEN 
PROBLEM 

The aim of this conclusion is twofold. We will first 
jive an overview of the different results presented in this 
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paper. Then we will discuss the problem of the size se- 
lection and the formation of barchans corridors. 

The starting point of the present work is the observa- 
tion that barchan dunes are organized in fields stretched 
along the dominant wind direction (figure 0). These 
barchans corridors are quite homogeneous in size and in 
spacing. For instance, the barchans field between Tarfaya 
and Laayoune presents the same five coherent corridors 
over at least 100 km. This size selection is of course 
not to be taken in a strict sense: there are large fluctu- 
ations from one dune to another, which have also to be 
explained. 

We have shown that the stability of a solitary dune 
essentially depends on the relationship between the size 
of the horns and that of the dune. Indeed, the dune 
receives at its back a sand flux proportional to its size 
but releases sand only by its horns. If the size of the 
horns is proportionally smaller for large dunes than for 
small ones, the steady state of the dune is unstable: it 
either grows or decay ( figure ITU|l . If, on the contrary, the 
sand leak increases faster than the dune size, it pulls the 
dune back to equilibrium. Furthermore we have shown 
that the fact that a dune is fed by the output flux of 
the dunes upwind does not change the stability analysis. 
This is essentially because a dune can influence another 
dune downwind through the flux but there is no feedback 
mechanism. 

We are thus left with a secondary question: how does 
the horns width evolve with the dune size? The only field 
measurements from which the horns size can be extracted 
are the shape measurements of eight dunes by Sauermann 
et al. [2]} . The sum of the width of the two horns is found 
to be between 12 and 28 to for the five small dunes they 
measured (2 and 3 to high), and between 12 and 17 to 
for the three larger ones (heights between 6 and 8 to). 
If we trust the relevance of their selection of dunes, this 
means that the horn size is almost independent of that of 
the dune. This is also coherent with their claim that the 
slip face is proportionally smaller and the horns larger 
for small dunes than for large ones. In that case, solitary 
barchans should be individually unstable. 

The second indication is provided by the C£ modeling, 
with which we recover that this steady state is in fact 
unstable (figures and 0) , in fact for the very same rea- 
sons as above. The solution can be artificially stabilized 
by putting at the back of the dune exactly what it looses 
by its horns but this is only a numerical trick. What 
determines the size of the horns in the model? The 3D 
solutions can be thought of coupled 2D solutions [28| . 
Then, the horns start when there is no slip face, i.e. 
when the length becomes of the order of the minimal 
size of dunes. This simple argument leads to think that 
the horns should keep a characteristic size of order of few 
saturation lengths L whatever the dune size is. 

We have shown that there is a second robust mecha- 
nism of instability. We know from field measurements, 
numerical models and theoretical analysis that the dune 
velocity is a decreasing function of its size. The reason 




FIG. 13: Aerial photographs of several barchan dunes in the 
same region as the field of figure Q They all exhibit an in- 
stability on the left side, leading to a periodic array of small 
slip-faces. 



is simply that the flux at the crest is almost independent 
of the dune size and will make a small dune propagate 
faster than a large one. This is sufficient to predict the 
coarsening of a dune field: because they go faster, small 
dunes tend to collide the large ones making them larger 
and slower. . . This collision instability should also lead to 
an ever growing big dune. 

The scales over which all instabilities develop are the 
relaxation time r r and length lambda r . For the eastern 
corridor of figure ^ t ne order of magnitude of the dune 
width is 100 to which gives r r ~ 1.1 century and A r ~ 
4 km. For the western corridors, the dunes are smaller 
(to ~ 20 m and the characteristic scales become r r ~ 
10 months and A r ~ 160 to. These lengthes are much 
smaller than the extension of the corridor (300 km) so 
that these instabilities have sufficient space to develop. 

However, the actual barchan corridors are homogenous 
and one cannot see any evidence of such instabilities. As 
a conclusion, the dune size selection and the formation of 
barchan corridors are still open problems in the present 
state of the art. There should exist another robust dy- 
namical mechanism leading to an extra-leak of the dunes, 
to balance the collisional and the flux instabilities. There 
are already two serious candidates for this mechanism. 
First, we do not have any information on the collision 
process. In the C£ model presented here, the collision 
of two dunes leads to a merging into a larger dune, but 
in the reality, it could lead to the formation of several 
dunes. Second, we have only investigated here the case 
of a permanent wind. We have shown that the dunes 
characteristic times are larger than one year so that the 
annual variations of the wind regime could have dras- 
tic effects. Figure IT51 shows aerial photographs of eight 
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barchan dunes. They all present an instability on their 
left side, leading to the formation of a periodic array of 
1 m high slip faces. This asymmetry suggests that it 
is due to a secondary wind (probably a storm) coming 
from the west (from the left on the figure). This insta- 
bility could lead to a larger time averaged output flux 
than expected for a permanent wind. Further work in 
that direction will perhaps shed light on the formation of 
nearly homogeneous corridors of barchan dunes. 



The authors wish to thanks S. Bohn, L. Quartier, 
B. Kabbachi and Y. Couder for many stimulating dis- 
cussions. 



Appendix: the 3D C£ model 

The three starting equations of the model are the con- 
servation of matter, the charge equation and the coupling 
between the saturated flux and the dune shape h: 

d t h + d x q = DV ■ (qVh), (32) 
OxQ = — ^ — , (33) 

^ - 1 + A [ ( dX . d x h e + Bd x h e . (34) 

We recall that the overall flux is the sum of q along the 
wind direction, plus an extra flux due to reptons along the 
steepest slope. The two last equations do not contain any 
y dependence and can thus be solved for each slice in x 
independently, using a discrete scheme in space (dx). The 
conservation of matter Ip?^ couples the slices through the 
diffusion term and is solved by a semi-implicit scheme of 
time step dt. To speed up the numerical computation of 
the saturated flux, we use the discrete Fourier transform 
T of the dune envelope h e : 

q S at = Q(l + T~ x {T{h e ){A\k\ + iBk)}) . (35) 

This envelope is composed of the dune profile h(x) up to 
the point where the turbulent boundary layer separates: 

x < xj, : h e (x) — h(x). (36) 

In the absence of any systematic and precise studies on 
this separation bubble, we assume that the separation 
occurs when the slope is locally steeper than a critical 



value fj,b = 0.25: 

h(xb) — h(x(, + dx) > fibdx. (37) 

When the dune presents a slip face, the boundary layer 
thus separates at the crest. The separation streamline is 
modeled as a third order polynomial: 

Xb < x < x r : h e (x) = a + bx + ex 2 + dx 3 . (38) 
The four coefficients are determined by smooth matching 
conditions: 

h e {xb) — h{xb), h e (xb — dx) = h(xb — dx), (39) 
h e (x r ) = h(x r ), h e (x r + dx) = h(x r + dx), (40) 

and the reattachment point x r is the first mesh point 
for which the slope is nowhere steeper than \Xb- There is 
no grain motion inside the recirculation bubble, so that 
the charge equation should be modified to d x q = —q for 
Xb < x < x r . Similarly, on the solid ground (h = 0) no 
erosion takes place, so that d x q = 0. 

The last important mechanism is the relaxation of 
slopes steeper than by avalanching. Rather than a 
complete and precise description of avalanches of grains, 
we treat them as an extra flux along the steepest slope: 

d t h + d x q = V- [(Dq + E5fi)Vh], (41) 

where Sfi is nul when the slope is lower than /i^ and equal 
to 6fj, — \Vh\ 2 — /it? otherwise. For a sufficiently large co- 
efficient E, the result of this trick is to relax the slope 
to /irf, independently of E. Note that, as the diffusion of 
reptons, these avalanches couple the different 2D slices. 
The value of the parameters have been chosen to repro- 
duce the morphological aspect ratios and are given in the 
following table: 



A = 


4.7 




curvature effect 


B = 


5.0 




slope effect 


D = 


0.1 




Lateral diffusion 


Mb = 


: 0.25 




separation slope 


E = 


50 




avalanches 


Hd = 


= 0.5 




avalanche slope 


dx = 


= 0.25 


- 1.0 


grid step 


dt = 


0.001 


— 0.1 time step 


M = 


: 64- 


512 


box size 



The results presented in this paper have been obtained 
for different discretization time and space steps, different 
box sizes and different total times. This explains the 
slight dispersion of the measurements. 
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